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A general analytical framework for the determination of the mean field states at arbitrary rational 
filling factors for the 2DEG in FQHE regime is given. Its use allows to obtain analytic expressions 
for the solutions at filling factors of the form u = 1/q for arbitrary odd q. The analysis can be 
performed for two general classes of states characterized by 7 = 1 or 7 = i particles per unit cell. 
Instead of the periodic peaks of the Wigner solid solution, the new states show electron densities 
forming percolating ridges that may favor an energy decrease through correlated ring of exchange 
contributions. Therefore, we estimate that they can realize mean field versions of the so called 
Hall Crystal (HC) states. The obtained analytic HC solution shows the same crystalline symmetry 
that the corresponding WC state in its class 7 = 1, but a qualitatively different charge density 
distribution. The energy dependence of the corresponding HC and WC states on the filling factor 
is also evaluated here for the class 7 = 1/2. The results show a crossover between HC state and the 
Wigner crystal, close to filling 1/7. Therefore, transitions may occur from one to the other as the 
electron density is varied. This result is consistent with recent experimental findings. 

PACS numbers: 73.43.Cd, 73.43.-f 



I. INTRODUCTION 

The quantum Hall effect has been the source of continued interest for already a quarter of a century,^ As the 
samples have become cleaner and the temperature made lower unexpected structures have been seen, suggesting a 
richer physics than originally thought. Although the system was believed to simply transit from a liquid state to a 
Wigner Crystal as the electron density is decreased, recent observations show that new structures develop that may 
signal transitions between states of different symmetryi^i^iii^ Numerical calculations also find a variety of states among 
which a system may transit, notably at Landau levels higher than the first i^ii^ 

The formation of charge density wave states in the thermodynamic limit may be studied using the Hartree-Fock 
approximation.^' The formalism allows for the consideration of states of various symmetries, unidirectional waves, 
triangular or square lattices, etc. The criterion normally used to choose among them is to select the state of lowest 
energy. Yet none of such approximate states is the true ground state, so the effect of correlations is a consideration to 
have in mind. In particular, depending on the distribution and relative phases of the occupied single particle states 
the correlations may be enhanced or suppressed. Since experiments show cases in which two different phases appear 
it is worthwhile to explore the various symmetries allowing for self-consistent solutions. 

In this work we generalize the formalism given in rcfcrence^"^ for the particular :/=l/3 filling, to perform an analytic 
treatment of the Hartree-Fock problem in the lowest Landau level (LLL) at arbitrary filling. We also will discuss in 
more detail particular solutions at filling factors v = 1/q for arbitrary odd q. For each value of q two classes of states 
are considered and indexed by the number of particles per unit cell 7=1 and 7 — 1/2. For each one of these two 
groups, two types of mean field states are associated. One in which the electrons build up well localized gaussian- 
like peaks: the Wigner Crystal (WC),and another in which the charge density forms hexagonal ridges increasingly 
resembling a bee hive as the electron density is decreased: the Hall crystal(HC) state. Therefore, the low density 
limit of this new state differs qualitatively from the one in the Wigner Crystal state a^^'^^i^^ . One can then expect 
that correlations arising from ring exchange will affect more this state, since percolating paths are open throughout 
the structure fi^ii^. For the particular 7 = 1/2 class, we find a crossover between the WC and HC state, showing that 
while the Wigner Crystal (WC) in the class has lower energy at fillings below 1/7, the new HC state however, has 
lower energy at larger fillings. Thus, even within mean field theory both states are competing as precursors to the 
ground state. 

Before ending this section, let us comment on the physical relevance of the parameter 7, the number of electron 
per unit cell. Its importance is related in one sense, with the fact that it determines the symmetry of the mean filed 
solutions. The value 7 = 1 is associated to the Yoshioka-Lee(YL) WC state which in the Hartree-Fock framework 
shows the lowest energy pere particle. Moreover, within the same class 7 = 1 there is a HC state which shows a bee 
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hive like density, and although having higher mean field energy than the YL state, one can suspect that it can receive 
highly cohesive correlation energy contributions. Further, the class 7 = 1/2, is the next in the, up to now, identified 
relevance, since for these states, it was made clear in Refj^, that although showing higher mean field energies than the 
YL states, its states predict energy gasps for all the relevant filling factors in FQHE. Since the HC in this class is also 
showing a bee hive like density distribution, the possibility that the correlation energies can make them competitive, 
a least as alternative ground states in the experimentally detected phase transitions, is a real one. a The possible 
relevance of mean field states having higher than two particle per unit cell is not excluded, but their consideration 
should wait for obtaining their wave functions. In this sense, we think that precisely the method outlined in this paper, 
can be generalized for obtaining 7 < 1/2 mean field solutions or at least a large amount of analytical information on 
them. 

In Sec. 2 a formula for the single particle Pock operator is presented for the case of an hexagonal lattice. Sec. 3 
is devoted to construct a basis in terms of which the Fock operator matrix representation reduces to q dimensional 
blocks, where q is the denominator of the filling fraction. Section 4 presents the exact determination of the new HC 
mean field state for each of the classes 7 = 1, 1/2. They are characterized by having a spacial region of very low 
charge density around a high symmetry point in the unit cell, that acquires a zero if the filling is of the form 1/q. The 
case q = 3 is discussed in more detail. The evaluation of the energy dependence on the filling factor for the WC and 
HC states showing a crossover, for the particular class of states 7 = 1/2 , is also done in this section. Finally Sec. 5 
presents our conclusions. 



II. THE FOCK OPERATOR 



We consider N^. electrons on a plane, in a strong perpendicular magnetic field. For such a system it was proven 
long ago that at filling less than one the only consistent mean field solutions have space fluctuations^^. We shall 
here assume the charge density to form a periodic lattice with a unit cell containing 7 electrons, this number being 
fractional or integer. Denoting by (j) the flux through such cell in units of the flux quantum (po = hc/e, one can readily 
verify the simple relation 

0j/ = 7. (1) 
The Fock operator in the LLL may then be written in the formic 



-ffffF = ^w(Q)exp(-^^^2^)Tr2„xQ, (2) 
Q 

where 

. .„„V(Q)e..(^2!) - ^y^)) jSl, (3, 

Here rg ~ ^ hc/eB is the magnetic length, Io{u) a modified Bessel function and Eq the background dielectric 
constant. The operator Ta displaces the function it acts upon in —a adding a magnetic phase factor, as defined in 
Appendix A. The Fourier components of the periodic charge density p(x) are defined as usual. 



P(Q) = ^ / dx p(x) exp(zQ.x), (4) 

Acell J 

where Aceii is the unit cell area 

Acell = n.aixa2 = 27rr^(;i, (5) 

n being a unit vector normal to the plane containing the electrons. Assuming triangular symmetry, the lattice formed 
by the electrons is then invariant under translations in the set of vectors 

R = niai+n2a2, ni, ^2 = 0, ±1, ±2, (6) 
ai=a(l,0), (7) 



A V3, /47r0 
a2 = a (-, — ), a = J -^fo- (8) 
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The vectors Q in Eq. (4)span all points in the reciprocal lattice and are given by 
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The Fourier components of the density obey the sum rul&i^ 

^ |2^r2p(Q)|2 exp(^) = - v). (10) 
Q ^ 

where the term Q = is omitted from the sum. This relation states that within the mean field approximation the 
liquid state of uniform density (p(Q) = all finite Q) in the LLL is only possible at filling one. At fractional fillings 
the right hand side is finite and so must be at least one finite wave- vector Fourier component of the charge density. 



III. BLOCK DIAGONALIZATION OF THE FOCK OPERATOR 



The mean field hamiltonian ([2]) describes an electron in a periodic potential and a perpendicular magnetic field, a 
case for which many results are knownJi An important property is that the single particle spectrum in the LLL is 
arranged in non overlapping bands, each with the same number of states.^? The number of bands equals the numerator 
of the flux per cell, the latter assumed a rational. For filling v = p/g, with p, q prime to each other, Eq.(IT]) yields a 
flux = 7g/p per plaquette. Assuming 7 to be a rational this flux is then also a rational number. 

We consider in what follows two simple cases that illustrate how different values of the flux ^ are to be treated. 
One is when this number is an integer, and another when it is half an integer. For simplicity we set 7 = 1, 1/2 and 
V — \/q, so that 4> = 9)Q/2, respectively. Our results cover the more general case 7 = p,p/2 and v = p/q, giving 
rise to the same values of the flux we include in the following discussion. Other cases may be treated using similar 
methods to the ones described below. 



A. Integer flux quanta per unit cell 

This case is important since it corresponds to the standard WC, in which each unit cell captures a full electron 
charge. The flux traversing a plaquette is q so that the single electron spectrum will have q bands, one of which is 
completely filled and the others empty. Since these bands do not overlapi^ the WC state has thus a gap for all values 
of q, whether even or odd. Because an essential feature of experiment at not too low filling fractions is the different 
behavior at even and odd values of such quantity, the WC state is not a good candidate for being the mean field 
precursor to the true ground state. 

Owing to definition ^ the magnetic translations entering the Fock operator have the form 

Tr^nxQ = T_Q2 ,Qi ■ (11) 

Since the flux piercing the unit cell is an integral number q of flux quanta, the set of translation operators Tr for all 
the R defined in ([5]) commute among themselves, allowing to find common eigenfunctions to all of them. This is not 
the case for the translations pT|) since the original unit cell is partitioned in smaller sectors if q is greater than 1. The 
basis we shall construct defines a set of g-dimensional subspaces, which are closed under the action of translations 
pT|) for all values of Q. 

A first step in finding the basis is to define a set of eigenfunctions Xk(x) of a translation in the vector —ai/q for 
each value of the momentum p=fik. Expressed as linear combinations of the functions (/?k(x) defined in Appendix A, 
we write them in the form 

2 

Xk(x)= C,(k)r_|ai¥'k(x), (12) 
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where, for definiteness, we have assumed q to be odd. These functions must obey the condition 

T_i Xk{x) = Axk(x). 

q ^ 

One finds for the eigenvalues A and coefficients Cs the set of q solutions 

AW(k)=exp(z^^+z — ) (13) 

, 1 , , sk.ai .27rrs 
k = — exp -I I 

2 2 ■ 

Substituting in (fT2)) yields the 5 eigenfunctions 

Xr(x) = — 5^ exp(-z— ^-*^)T_|a, ¥'k(x), (14) 

^~ 2 

k = k + nsi+ms2. n, m = 0, ±1, ±2, ... (15) 

The last relation expresses the fact that the states in the new basis are equivalent upon a shift of k in any linear 
combination with integer coefficients, of the unit cell vectors of the reciprocal lattice corresponding to the periodicity 
of the density. The equivalence follows from the following properties: (a) the functions lyfk(x) are eigenfunctions of 
any translation Tr for lattice vectors R given by ([5]) , (b) the operator Tr commutes with all translations entering in 

the definition of x^'^-'(x), and (c) relation (jA8|) in Appendix A, stating the equivalence between magnetic translations 
acting on <Pk(x) and a shift in the momentum labelling these functions. 

Let us now inspect the effect of a magnetic translation in ai.2/ q on the new fimctions. If such a transformation leaves 
the g-plets invariant, then the matrix reduction of the Hartree-Fock hamiltonian will follow. One has, 

T, J^\^) = — J2 exp(-z— ^-z--)n,^r_..,^k(x). (16) 

^~ 2 

After using (jA4p for changing the order of the two operators within the sum, it follows that 

T,^J^\^) = exp(-z il^)x(['^-^l'(x), (17) 

where the square bracket defines the number in the set { — -^^^2^1 '^^^^I which is equivalent, modulo q, to the integer 
in the argument. Thus, a magnetic translation in ^0.2 just turns one function in the Q'-plet into another. Besides the 
properties already discussed, the basis can be checked to obey 



7aiXr(x) = exp(-z k.ai )xr(x), (18) 
T.,x'^\^) = exp(-* k.a2 )Xk''^(x), (19) 
Pxi'\^)^xi-^\^), (20) 



where the parity transformation P is defined as usual, P Xk^''(x) = Xk^''(~x). 

^From the above considerations it follows that the q-dimensional subspace spanned by the functions Xk^"* (x) at fixed 
values of k is left invariant by the action of the operators pT|) for arbitrary values of Qi and Q2- Since the Fock 
hamiltonian involves just a sum of such translations, it leaves invariant these g-dimensional subspaces, as well. It is of 
interest to note that this is a purely kinematic result which does not depend on the form of the interaction potential. 

In the new basis, the q^ matrix elements of the hamiltonian ([2|) can be readily found to have the convenient form 



(r',' 



<xi''^\HHF\x'; 



(r) 



E 



> 

w(Q)exp( — ^)exp 



(21) 



-ik.n X Q rl + i-Q2{2r 
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The problem has thus been reduced to the self-consistent diagonalization of a gr-dimensional matrix for each value of 
the wave vector k. For a sample of surface S the degeneracy D — BS/4>o of the Landau level of the non interacting 
problem is then split into q bands that span their range as k covers the Brillouin zone, each holding exactly D/q 
single-particle states. 



B. The 7 = 1/2 class 

We next turn our attention to the case v — l/q,j ~ 1/2 for which, following Eq. ([T|), cf) = q/2. This case is 
particularly interesting because if q is even, say q = 2r, then (j) — r and the single particle spectrum has just r = q/2 
bands. Since v = (l/2)/r only the lowest of these energy bands has occupied states, yet is only half filled. The 
Fermi energy is at the center of the band and the state is metallic. By contrast, if q is odd there are q bands in 
the spectrum, one of which is completely filled and the others empty, leaving the Fermi level in a gap. Even and 
odd filling fraction denominators thus show qualitatively different behavior, a metal or an insulator, as experiment 
demands. This remarkable property makes this state a reasonable candidate to be the mean field precursor to the 
true ground state of the system. 

The magnetic translations defining the Fock operator have now the form, 

Trl nxQ = r_|Q.^^^_^|Q^^^. (22) 

We define a doublet invariant under magnetic translations in the vectors ai and a2, 

(^^(x) = ^((Pk(x) + -— ^-—Ta,(pk(x)), a = ±1. (23) 
V2 exp(-ia2.k) 

These functions obey 

Ta,(/p^(x) = exp(-zai.k)(^-'^(x), (24) 
T'as'PkW = crexp(-ia2.k)(/3k(x), 

being turned into each other by translations along the axes, save for a phase factor. In analogy with the previous 
subsection we define next a set of eigenf unctions of the translations Ta , with the form 

xr^(x) = ^ ± exp(^^-.^)^a..^^£(x), (25) 
y/q q q ' 

^ — 2 

k = k + rt Si/2 + m S2/2, n, to = 0, ±1, ±2, ... 

Note that the last line indicates for this case that the equivalence of states is now under shifts in half the reciprocal 
lattice unit cell vectors. This is related to fact that magnetic translations in vectors (6) are non commuting, so that 
they have no common eigenfunctions. However, translations in twice the spatial unit cell vectors, are commuting 
operations. Similarly as in the previous subsection, the equivalence follows after considering that, (a) the functions 
(/?k(x) are eigenfunctions of any operator T^r. for lattice vectors R given by (b) T2R commutes with all the 
translations entering in the definition of Xk^''^'(x) through the original functions </3k(x) and (c) relation (|A8p in 
Appendix A, expressing the equivalence between magnetic translation on (pk(x), with a shift in the momentum 
labelling these functions. 

Aside from relations (|24p these functions satisfy the following transformations 

r|a.xi:"'(^)=exp(z^+z^)x^^\x), (26) 

r|a.xi:"^(^) =exp(-^^k.a2)xL'^-''"'(x). (27) 

Consider now fixed values of the quantum numbers (k, a). Relations (1^5)1 and P7)) directly show that all translations 
(1^^ included in the Fock operator ^ leave invariant the q-dimensional subspace spanned by the set X]^'''\x), all r. 

Using the commutation properties (jA4p and relations ([^5)1 and ([77)1 the matrix elements of the hamiltonian ^ can 
now be written as 



6 



^2q2 

Q 



„ 9 lnQi , 
-ik.n X Q + i — —{r + i 



l)^ ^r',[r+2Qi] 



(28) 



Again the mean field problem has been reduced to the diagonalization of a g-dimensional matrix, but now this must 
be done for each value of the wave vector k and the index a. An interesting outcome is that the matrix representing 
the hamiltonian is identical for the two values of cr, so that its q eigenvalues are twice degenerate. The obtained block 
diagonalization of the problem directly indicates a procedure for obtaining approximate selfconsistent solutions. One 
starts assumm=ing an initial density function showing the periodicity in the lattice R and perform its Fourier series 
in terms of the momenta in the reciprocal. This furnish the initial Fourier component of the density p(q) lattice cell. 
Then, defining a 2D uniformly spaces partition of the sets of vectors q for which the HF states are inequivalent, the 
set of Hamiltonian matrices associated to these momenta values are diagonalized. Note that these matrices are fully 
defined after specifying the initial Fourier component of the density. Then, the information in this diagonalization 
process furnish the values of the eigenfunctions coefficients. Next, the density to be employed in a new step can be 
calculated from the formula 



('?-l)/2 
cr=±l r=-((;-l)/2 



(29) 



The coefficients gr(k) are the lowest energy eigenvectors of the reduced Hamiltonian block matrix for each k value. 
Since these matrix blocks do not dependent of a the eigenvectors are .9r(k), r — {q — l)/2^ .... — {q — l)/2 are also a 
independent. That is, the same is to be used for both values of a in the 7 = 1/2 case. Further, this new density 
value is again Fourier expanded in the reciprocal lattice vectors and the Fourier coefficients employed to define new 
Hamiltonian block matrices to be diagonalized. Then, the described process should be developed iteratively up to the 
arrival to a convergence. To a situation in which the density at any new steps closely approaches the one the previous 
step. 



IV. ANALYTIC FORM OF THE HC STATES AT v = l/q FOR THE TWO 7 = 1, f CLASSES 

As discussed at the end of the last section, further progress will normally require a numerical routine that diago- 
nalizes self-consistently cither (PT|) or (1^5)) . where the solutions obtained at the end generate the same potential w(Q) 
that gave rise to them. There is a special situations, however, for which further analytic progress is possible. Because 
these sates qualitatively differs from the one in the usual WC, in a way suggesting the connection with the states 
discussed in^^ we call it below the Hall Crystal (HC) state. 

Its characterization rests on the following observation. Earlier numerical work showed that in the special case 
V — Xjq there are two pairs of special self-consistent solutions, one pair for each class of states 7=1 or 7 = 1/2 
respectiveljJ^. For a given class (7 = 1 or 7 = i) one of the states in its corresponding pair represents the WC. This 
wave-function shows gaussian-like peaks of its density centered at equivalent symmetry points in the lattice. However, 
in the second state in the pair, named before as the Hall Crystal (HC), the charge density vanishes at a single point 
in the unit cell and all its periodic images in the rest of the crystal. Further inspection in the vicinity of such points 
showed this zero to be of order 27(5 — 1). Since the total particle density is a sum of the positive definite particle 
densities of the single particle occupied states, all associated orbitals must vanish at the special points as the power 
7(q — 1). This information was extracted form the numerical solutions of the problemsi^. Below, we will use this 
condition to analytically determine the form the particular HC states. 

In terms of the basis states defined previously the eigenfunctions of the Fock operator may be written in the form 



(30) 



(31) 
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where 6 = 0, 1, g — 1 is the band index, and the label o = ±1 is to be omitted if 7 
in the LLL, in our sign convention (e = — |e|) they must be of the formic 



1. Because these functions are 



vI/(z,z*) = F(z*)exp(-^), 

where z = x + iy, z* = x — iy, and F{z*) is an analytic function of its argument. With no loss of generality we 
choose one of the special zeroes to be at the origin. As one approaches this point one expects that asymptotically 
F{z*) ~ {z*y, with s = 7(g — 1). Thus the function ^ itself and its first s — 1 derivatives must vanish at the origin, 
giving in all s independent equations to be satisfied. These, together with the normalization condition (j3ip total s + 1 
equations, sufficient to determine the q coefficients {gj!(k),r = — of the occupied orbitals. Specifically, 
for 7=1 one has s + 1 = q, while for 7=1/2 one has 2s + 1 = q, the factor of two arising from the double valued 
index a. Thus the full structure of the single particle wave functions are determined after solving the considered set 
of equations for the coefficients. 

Let us examine in more detail the two states for the class 7 = 1 and the simpler case g = 3. The WC has been 
reported in the literature for this case,— but the HC has not. Solving for the coefficients 55^.]^(k), (7g(k), ^^(k) defining 
the filled orbitals and substituting them in the equation 



E 5j'(k)xL'-'(x 



(32) 



where the first sum runs over all momenta k in the Brillouin zone, the particle density may be obtained. The result is 
shown in Fig. f (b), where the density for the Wigner Crystal solution at the same filling. Fig 1 (a), has been included 
for comparison. Notice the presence of sharp hexagonal channels surrounding low density regions where the density 
vanishes as the fourth power of the distance. The fact that the density percolates the structure much like the wax in 
a bee hive marks the essential difference between the HC and WC states. In the latter the charge density is made up 
essentially of gaussian functions centered at lattice points as seen in Fig. 1(a). As the filling fraction decreases these 
gaussian peaks become sharper and the limit of a classical point-electron WC is approached. 




(a) 



(b) 



FIG. 1: Particle density for the (a) WC and (b) HC states at filling = 1/3 and 7 = 1 (one electron per cell). Notice the 
percolating hexagonal ridges and the absence of any gaussian-fike peaks in the latter case. 



Further, insertion of the calculated density in Eq. ^ and use of this result in Eq.® allows finding the associated 
eigenvalues e'', r = 0, 1, 2 by diagonalizing the 3x3 matrix ((2T|) . Three bands are obtained, that span their range as k 
covers the Brillouin zone. 

The bands dispersion relations are illustrated in Fig. 2. They are quite narrow, with the lowest (the filled one) well 
separated from the rest by a sizable gap. The same pattern was found for larger values of q. The LLL having been 
split into q separate bands appears to yield the spectrum associated with a value of the magnetic field reduced by 
a factor l/q. This feature is reminiscent of the composite fermion theory, which interprets the FQHE as the integer 
QHE of composite fermions in a field reduced by the same factori^i^i 

The energy per particle of the HC solution we have just discussed is found through 



, 1 V^ffl = _o.362A 

N ^ 2 Eoro 

k 



(33) 
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FIG. 2: Dispersion relation for the three bands into which the lowest Landau level is split by the action of the Coulomb 

2 

interaction, for the same case as in Fig. 1. Energy units are y-^- 

The value obtained is 7% above the WC solution associated to the 7=1 class, for which case e = — 0.388 j^^. The 
usual minimal energy criterion thus points to the WC as a better candidate for the ground state and the HC an 
exited state, a conclusion that persists for larger values of q. Still, how these only approximate mean field solutions 
are affected by correlations is largely unknown save for second order perturbation corrections. 




FIG. 3: Energy per particle for different values oi q = \lv for the Wigner Crystal (dashed line) and the Hall Crystal (broken 
line) both for the class 7 = 1/2. Units are as in Fig. 2. Lines are drawn to guide the eye. 



A more interesting case corresponds to the class 7 = 1/2 for q=3. A similar analytic procedure as employed above 
gives the energy for both the HC states. However, It should be noticed that this special solution was already discussed 
in referenceiS. Concretely, the set of three equations 

1 

^ .g::(k)x^^Ho)-0, a = ±l (34) 

r=-l 
1 

^ g:\^)gl{^) = 1, (35) 

r=-l 

for the coefficients of the functions defining the HC states through its zero should now be solved. Note that in this 
case 7=1/2 the three equations correspond to the vanishing of the functions ((30)) at the origin for the two values of 
cr=±l plus the normalization condition (|3ip . 

Finally, the dependence of the energy per particle on the values of q for the class 7=1/2 was evaluated, for the 
WC as well as for the HC states. The results are shown in Fig. 3. Notice that the curves intersect at about filling 
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1/7, suggesting that phase transitions may occur between the HC and the WC states associated to the class having 
half particle per unit ceU (7 — 1/2), as the magnetic field is varied. It should stressed that, although the mean field 
results for the energy in the class of states 7 = 1/2, depicted in Fig. 3 are higher than for the ones in the class 7 = 1 
( -0.388 e^/ro aX v — 1/3 for the WC state in this class), these states could yet enter to compete in having lower 
energies after the inclusion of the correlations. Moreover, the fact that recent experiments detect phase transitions in 
varying the filling factors as well as the temperature, could evidence these states as alternative excited configurations 
which could become the ground ones after the temperature and the filling factor are varied. These possibilities are 
consistent with the experimental data reported in Ref.^. The evaluations of the energy per particle, for the WC and 
HC states of class gamma=l/2, were done by numerically solving the HF problem. The numerical solutions were 
obtained by using the same programs employed in Ref J^, which implements an algorithm being equivalent to the one 
sketched at the end of Section 3. In fact the analytical results presented in Sections 2 and 3, can be interpreted as 
obtained in order to perform in an analytic way a great deal of the numerical steps in the programs employed for 
obtaining the results in Refi^. 

V. SUMMARY AND CONCLUSIONS 

We have shown that analytical solutions to the Hartree-Fock problem for a two dimensional electron gas at filling 
factors V = p/q may be obtained for charge density wave states whose unit cell include 7 electrons, with this number 
fractional or integral. Our results follow from symmetry considerations and the construction of a special set of common 
eigenfunctions to all magnetic translations in a lattice vector. The spectrum of single particle states is shown to be 
organized in multiplets of dimension q in such a basis, which are left invariant by the action of the Fock operator. 
The matrix elements representing the Fock operator in each multiplet are explicitly determined. 

We treat in detail the case v = 1/g for the two classes holding one (7 = 1) or a half electron (7 = 1) per plaquette. 
In particular for each of such classes, we analytically determine a special state whose low density limit is not a Wigner 
Crystal. We call this state a Hall Crystal. It is characterized by occupied orbitals all of which have a zero of order 
"f{q— 1) at a point of high symmetry in each cell. The charge density is made of percolating hexagonal ridges expected 
to provide easy paths for correlated ring exchange that would lower the energyJ^ For the particular case, 7 — 1/2, the 
mean field energy was evaluated here as a function of the filling factor and shows a cross-over from a regime in which 
the Wigner Crystal like state has lower energy, to one in which the Hall Crystal one has lower energy. In contrast, 
the same evaluations for the class 7=1 down to fillings as low as 1/9, indicates that the Wigner Crystal like state 
has the lowest energy within our mean field approximation. The correlation energy affecting in a different way both 
states may alter this ordering, however. 

In concluding let us expose some perspectives for the extension of the work. Let us recall states in the lowest 
Landau level obeying periodic boundary conditions must vanish at a number of points equaling the degeneracy of 
the non-interacting systemii^ For v = 1/q this number is qN^, where iVg is the total number of electrons in the 
sample. But, for example, the Hall Crystal states at 7 = 1 attaches {q — l)Ne of them to fixed high symmetry periodic 
points in the electron lattice, leaving Nf^ zeroes whose spatial location depends on the quantum number k labelling 
each occupied state. Also, the associated determinant will vanish as the first power as two particles approach each 
other. However, the Laughlin state does so as a power q. This feature decreases the direct Coulomb energy in the 
Laughlin states by keeping particles as far as possible from each other, exhibiting built-in correlations which mean 
field solutions lack. However, the just described property opens the interesting idea of constructing new correlated 
states optimizing short as well as long range correlations starting form the HC states discussed in this paper. For 
this purpose, we suggest here to substitute the common factor of radial functions of the Nf. particles which appear 
in the HC state determinant, by Jastrow factors having the same number of zeroes that factor. In this way, the so 
constructed wave-functions will show q zeros when any two particle joins. But moreover, it can be expected that 
they will even more optimize the energy, since they also incorporate the long range crystalline correlations present in 
the Hall Crystal states. We hope the presented results and ideas for the extension of the work will stimulate further 
research directed to check whether the corrections to the mean field approach early introduced by one of the authors 
(F.C) could furnish a general theory for the FQHE. 

Acknowledgments 

We acknowledge partial support from Fondecyt, Grants 1060650 and 7020829, the Catholic University of Chile 
and the Caribbean Network on Quantum Mechanics, Particles and Fields of the ICTP Office of External Activities 
(OEA). One of us (F.C.) would like to thank the ICTP for hospitality while this work was being done. 



10 



APPENDIX A: EIGENFUNCTIONS OF MAGNETIC TRANSLATIONS 

It is well known that Bloch-like states may be constructed by placing a seed function at each lattice point and 
attaching to it an appropriate phase factor. In the presence of a magnetic field this procedure may be implemented 
using the zero angular momentum eigenfunction 

1 , x2 



and forming the sumi 




iVk = / >: (-1) ' ' exp(i kJ--2 ) 



The summation indices £i,£2 run over all integers, defining a planar lattice L through £ = ^ibi + t2^2 with the unit 
cell intercepting one flux quantum, so that n.bixb2 = ZTrr^. Displacements of the seed function are effected by 
magnetic translation operator Ta in the vector a, whose action on any function / is defined hy~ 

Ta /(x)=exp(-A(a).x) /(x-a). (A3) 

Here the vector potential is assumed in the axial gauge A(x) — B{~X2,Xi,0)/2 and the electron charge e is taken 
with its negative sign (e = — |e|). These translations in general do not commute, 

2ie 

T^,T^, = exp(— A(ai).a2)Ta,Ta„ (A4) 
= exp(— A(aJ.a2)rai+a2. 

In the special case of displacements in any vector belonging to L, however, since the flux trapped by any parallelogram 
bounded by lattice vectors is an integral number of flux quanta, all translations commute. 

One can easily check that the functions are eigenstates of translations in any lattice vector, satisfying the 
eigenvalue equation 

T, ^k(x) = Ak(^) (/5k (x), (A5) 
AkW = {-if'^^ exp(-i k.^). (A6) 

Arranged in a Slater determinant these functions are exact solutions of the Hartree-Fock proble m^^i^^i^'^i^^ . This 
strong property arises from the fact that the HF single particle hamiltonian commutes with all translations leaving L 
invariant. The functions (|A2[1 are common eigenfunctions of the commuting magnetic translations. Moreover, the set 
of eigenvalues (|A6[) uniquely determines them. Therefore, the HF hamiltonian associated with the Slater determinant 
can not change those eigenvalues and the <y9k should be eigenfunctions. 

An important property of the basis functions (|A2[) is that an arbitrary translation is equivalent to a shift in the 
momentum label, modulo a phase factor. Operating twice with the translation operator involving an arbitrary vector 
a and a vector in the lattice t, and using Eqs. (|A4|1 and (|A5|) one readily gets, 

TaT, ^k(x) = AkWra^k(x), (A7) 

= exp(^A(a).^) T.Ta (^k(x), 
nc 

which can also be written as 

TfTa (^k(x) = Ak+2|A(a)Wra <^k(x). 

Then, taking into account that the set of eigenvalues defines uniquely the wave-functions modulo a phase, it follows 
that 

Ta ¥'k(x) - .Fk(a) ^k+||A(a)(x), (A8) 
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